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Abstract. We have performed evolutionary calculations 
of very-low-mass stars from 0.08 to 0.8 Mq for different 
metallicites from [M/H]= -2.0 to -1.0 and we have tabu- 
lated the mechanical, thermal and photometric character- 
istics of these models. The calculations include the most 
recent interior physics and improved non-grey atmosphere 
models. The models reproduce the entire main sequences 
of the globular clusters observed with the Hubble Space 
Telescope over the afore-mentioned range of metallicity. 
Comparisons are made in the WFPC2 Flight system in- 
cluding the F555, F606 and F814 filters, and in the stan- 
dard Johnson-Cousins system. We examine the effects of 
different physical parameters, mixing- length, a-enriched 
elements, helium fraction, as well as the accuracy of the 
photometric transformations of the HST data into stan- 
dard systems. We derive mass-effective temperature and 
mass-magnitude relationships and we compare the results 
with the ones obtained with different grey-like approxima- 
tions. These latter are shown to yield inaccurate relations, 
in particular near the hydrogen-burning limit. We derive 
new hydrogen-burning minimum masses, and the corre- 
sponding absolute magnitudes, for the different metallici- 
ties. 

We predict color-magnitude diagrams in the infrared 
NICMOS filters, to be used for the next generation of 
the HST observations, providing mass-magnitudes rela- 
tionships in these colors down to the brown-dwarf limit. 
We show that the expected signature of the stellar to 
substellar transition in color-magnitude diagrams is a se- 
vere blueshift in the infrared colors, due to the increasing 
collision-induced absorption of molecular hydrogen with 
increasing density and decreasing temperature. 

At last, we apply these calculations to the observed 
halo field stars, which yields a precise determination of 
their metallicity, and thus of their galactic origin. We find 
no evidence for significant differences between the halo 
field stars and the globular cluster sequences. 



Subject headings : stars: Low-mass, brown dwarfs 
stars: evolution — stars: globular clusters 



1. Introduction 

Over the past decade considerable effort, both observa- 
tional and theoretical, has been directed towards a more 
accurate determination of the stellar lower main sequence, 
down to the edge of the sub-stellar domain. Such a de- 
termination bears important consequences for our under- 
standing of a wide variety of astrophysical problems, from 
star formation and stellar structure to galactic formation 
and evolution. Although the lower main sequence of the 
disk is relatively well determined since the survey of Monet 
et al. (1992), the situation is less well defined for the Galac- 
tic halo, primarily because of the greater difficulties in- 
volved in identifying the halo stars. The task of detecting 
halo low-mass stars (LMS) and measuring their magni- 
tude is of formidable difficulty with ground-based tele- 
scopes. Globular clusters (GCs) have always presented a 
particular interest for the study of the stellar halo since 
we can more easily determine their main sequence than 
for field halo stars. Most of the observations of GCs fo- 
cussed on the upper main sequence, i.e. the turn-off point, 
and the red giant branch, for a comparison of this region 
with theoretical isochrones yields a determination of the 
age of the clusters, and thus a lower-bound for the age of 
the Galatic halo. Thanks to the tremendous progress in 
deep photometry realized recently with the Hubble Space 
Telescope (HST), which reaches unprecedented magnitude 
and spatial resolution, the lower main sequence of GCs is 
now observed nearly down to the hydrogen burning limit. 
Thanks to the high angular resolution achievable with the 
HST, accurate photometry is feasible to levels about 4 
magnitudes fainter than with ground-based observations, 
allowing photometry of very faint stars. Several HST ob- 
servations of globular clusters are now available, spanning 
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a large metallicity-range from solar value to substantially 

metal-depleted abundances, as will be presented in the 
next section. The lower main sequence of these clusters 
is well defined and offers a unique possibility to probe 
low-mass star evolutionary models for various metallici- 
ties down to the hydrogen-burning limit. 

In spite of considerable progress in stellar theory - in- 
ternal structure, model atmospheres and evolution - all 
the LMS models so far failed to reproduce accurately the 
observed color-magnitude diagrams (CMD) of disk or halo 
stars below ^ 4000 K, i.e. ^ 0.4 — 0.6 A^q, depending on 
the metallicity. All the models predicted too hot an effec- 
tive temperature for a given luminosity, i.e. were too blue 
compared to the observations by at least one magnitude 
(see e.g. Monet, 1992). Such a disagreement stemed essen- 
tially from shortcomings both in the physics of the inte- 
rior, i.e. equation of state (EOS) and thus mass-radius re- 
lationship and adiabatic gradient, and in the atmosphere, 
since all models were based on grey atmospheres and a dif- 
fusion approximation. Important progress has been made 
recently in this field with the derivation of an appropri- 
ate EOS for low-mass stars and brown dwarfs (Saumon, 
Chabrier & Van Horn 1995; see Chabricr and Baraffe 
1997), non-grey model atmospheres for M-dwarfs (AUard 
and Hauschildt 1995, 1997; Brett 1995) and evolutionary 
models based on a consistent treatment between the in- 
terior and the atmosphere profile (Baraffe, Chabrier, Al- 
lard & Hauschildt 1995; Chabrier, Baraffe & Plez 1996; 
Chabrier & Baraffe, 1997). These models now reach quan- 
titative agreement with the observations, as shown in the 
afore-mentioned papers and presented in §4. 

Globular clusters offer the great advantage of all stars 
having the same metallicity, determined relatively accu- 
rately from bright star spectroscopic measurements. More- 
over they are old enough for all the stars to have reached 
thermal equilibrium, so that age effects do not affect the 
luminosity of the objects near the bottom of the MS. For 
these reasons, the mass-luminosity relationship is well de- 
termined along the entire MS of GCs, from the turn-off 
down to the brown dwarf limit, with no dispersion due 
cither to age or metallicity. From the theoretical view- 
point, these properties restrain appreciably the degrees of 
freedom in the parameter space, so that globular clusters 
provide a very stringent test to probe the validity of low- 
mass star evolutionary models and of the related mass- 
luminosity relationships for various metallicities. Agree- 
ment between these models and observations is a neces- 
sary condition (but not sufficient !) to assess the validity 
of these mass-luminosity relationships, a cornerstone to 
derive reliable mass-functions. 

In this paper, we present new evolutionary models 
for metal-depleted {[M/H] < —1) low-mass stars (m < 
0.8 Mq), based on the most recent non-grey model atmo- 
spheres. We compare the results with the observed CMDs 
of three globular clusters, namely NGC6397, M15 and 
u)Cen, for which HST observations are available. The pa- 



per is organized as follows : the observations are summa- 
rized in §2 whereas the theory is outlined in §3, where 
comparison is made with other recent LMS models. Com- 
parison between theory and observation is presented in §4. 
Section 5 is devoted to discussion and conclusion. 

The extension of the present calculations to more 
metal-rich clusters ([Fe/H] > -0.5) and field stars, which 
require more extensive calculations, will be presented in 
a forthcoming paper (AUard et al. 1997b), as well as the 
derivation of the mass function for globular star clusters 
and halo field stars, from the observed luminosity func- 
tions (Chabrier & Mera, 1997). 

2. Observations 

The decreasing luminosity of stars with mass {L ~ 
j^2.5-3.5 -f-Qj. low-mass stars) renders the observation of 
the lower main sequence almost impossible from ground- 
based telescopes. Despite these difficulties, a number of 
clusters have been investigated with large telescopes by 
several groups worldwide to determine the bottom of the 
main sequence. To our knowledge, the most extensive 
study was conducted by Richer et al. (1991) who deter- 
mined accurately the MS of six GCs up to a magnitude 
V ^ 24; My ~ 10, i.e. a mass m 0.4 M©. Although in- 
trinsically interesting, these observations do not provide 
any information about the very bottom of the MS, and 
the shape of the luminosity function near the hydrogen- 
burning limit. Moreover, for metal-depleted abundances, 
the afore-mentioned limit magnitude and mass correspond 
to a temperature Tes ^ 4000 K (see §4). Above this tem- 
perature, the physics of the stellar interiors and atmo- 
spheres is relatively well mastered, so that these observa- 
tions put little constraint on the models, and the related 
mass-luminosity relationship. 

A recent breakthrough in the field has been accom- 
plished with the deep photometry of several GCs obtained 
with the refurbished Space Telescope. Several observations 
are now available for NGC6397 (Paresce, DeMarchi & 
Romaniello, 1995; Cool, Piotto & King, 1996), NGC7078 
{= M15) (DeMarchi & Paresce 1995; Piotto, Cool & King, 
1996)and coCen (Elson, Gilmore, Santiago & Casertano, 
S., 1995). Observations for these clusters reach V ~ 27, 
/ ^ 24, almost the very bottom of the main sequence, and 
thus provide a unique challenge to probe the validity of 
the LMS models down to the brown dwarf regime. 

• Metallicity. All the afore-mentioned clusters arc sub- 
stantialy metal-depleted. Before going any further, it is es- 
sential to define what we call "metallicity" in the present 
context, for different definitions are used in the literature. 
What is observed in globular clusters is the iron to hydro- 
gen ratio [Fe/H]. The continuous production of oxygen 
in type II supernovac during the evolution of the Galaxy 
has led to the well observed enhancement of the oxygen to 
iron [O/Fe] abundance ratio in old metal-poor stars with 
respect to the young disk population. Since our basic mod- 
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els assume a solar-mix composition, i.e. [Fe/0] = 0, the 
afore-mentioned oxygen-enhancement must be taken into 
account to make consistent comparison between theory 
and observation. We use the prescription derived by Ryan 
and Norris (1991) for halo subdwarfs, i.e. 

[M/H] « [O/H] = [0/Fe] + [Fe/H] (1) 

with 

[0/Fe] = +0.35 for [Fe/H] < -1 

[O/Fe] = -0.35 X [Fe/H] for - 1 < [Fe/H] < 

Thus a cluster with an observed [Fe/H] = —1.3, 
for example, corresponds to a model with a metallic- 

ity [Z] = log(Z/Z0) = [M/H] 1.0, certainly not 

[M/H] = —1.3. This latter choice does nut take into ac- 
count the enhancement of the a— elemental and yields in- 
consistent comparisons, especially at the bottom of the 
MS where the stellar optical spectra and colors are shaped 
by H~, MgH, CaH and TiO opacities. The Ryan & Nor- 
ris correction is based on the spectroscopically-determined 
abundance of 370 kinematically selected halo stars in the 
solar neighborhood. Although this procedure is not as fully 
consistent as calculations conducted with the exact mix- 
ture, it certainly represents a fairly reasonable correction 
for the metal-poor star oxygen enrichment. The accuracy 
of this prescription will be demonstrated in §3.4, where 
calculations with the appropriate a-enhanced mixture are 
presented. 

• Photometric conversion. The afore mentioned clus- 
ters have been observed with the Wide Field and Plane- 
tary Camera-2 (WFPC2) of the HST, using either of the 
FilAW - /, F555M^ - F or F606M^ filters, where /, V 
refer to the standard (Johnson-Cousins) filters. Thanks to 
the courtesy of A. Cool, I. King, G. DeMarchi, F. Paresce, 
G. Gilmore and R. Elson, we have been able to use the 
data in these filters, in the so-called WFPC2 Flight system 
for A''GC6397, M15 and uiCen. For these clusters, com- 
parison between theory and observations is made directly 
in the Flight system, thus avoiding any possible uncertain- 
ties in the synthetic flight-to-ground photometric transfor- 
mations of Holtzman et al. (1995). The model magnitudes 
were calculaied using the filter transmissions curves and 
the observect^ zero points prescribed by Holtzman et al. 
(1995). 

The conversion of apparent m (observation) to abso- 
lute M (theory) magnitudes requires the knowledge of 
the distance modulus (m — M), corrected for the inter- 
stellar extinction in each filter. In order to minimize the 

* The a-elements include 0,Ne,Na,Mg,Al,Si,P,S,Cl,Ar,Ca 
and Ti 

** Since observed zero-points were not available for the F606W 
filter, we used the synthetic value listed in table 9 of Holtzman 
et al. (1995) for this filter. An inspection of this table confirms 
that synthetic and observed zero points agree almost exactly 
for the standard WFPC2 filters. 



bias in the comparison between theory and observation, 
we have used the analytical relationships of Cardelli et 
al. (1989) to calculate the extinctions from the M-dwarf 
synthetic spectra of AUard & Hauschildt (1997) over the 
whole frequency-range, from the reddening value E{B—V) 
quoted by the observers, and we have compared the ob- 
served data with the theoretical models corrected for red- 
dening in the WFPC2 filters, when available. The extinc- 
tion in each filter was found to depend very weakly on the 
spectral type (^ 0.05 mag), and thus to be fairly constant 
over the whole considered temperature range. We believe 
these determinations, based on accurate synthetic spec- 
tra, to yield the most accurate extinction and reddening 
corrections and the most consistent comparison between 
theory and observation. These values are given in Table 
I, for each cluster, and compared to the values quoted 
by the observers. Table I also summarizes the character- 
istics adopted for the three clusters of interest. Note that 
some undetermination remains in the distance modulus of 
NGC6397, yielding a difference in the magnitude of ^ 0.2 
mag, as will be shown in §4. This stresses the need for a 
more accurate determination of this parameter. 

3. Theory 

We have derived recently evolutionary models aimed at 
describing the mechanical and thermal properties of LMS. 
These models are based on the most recent physics char- 
acteristic of low-mass star interiors, equation of state 
(Saumon, Chabrier & VanHorn 1995), enhancement fac- 
tors of the nuclear rates (Chabrier 1997) and updated 
opacities (Iglesias & Rogers 1996; OPAL), the last gener- 
ation of non-grey atmophere models (AUard & Hauschildt 
1997, AH97) and accurate boundary conditions between 
the interior and the atmosphere profiles. This latter con- 
dition is crucial for an accurate description of LMS evo- 
lution, for which any diffusion approximation or grey 
treatment yield inaccurate results (Baraffe, Chabrier, Al- 
lard & Hauschildt 1995; Chabrier, Baraffe and Plez 1996; 
Chabrier & Baraffe, 1997). A complete description of the 
physics involved in these stellar models is presented in 
Chabrier and Baraffe (1997) and we refer the reader to 
this paper for detailed information. 

Below ~ 0.4 Mq, the stellar interior becomes fully con- 
vective : the evolution of LMS below this limit is rather 
insensitive to the mixing length parameter a = l/F[p, 
and thus the models are not subject to any adjustable 
parameter. From this point of view, VLMS represent a 
formidable challenge for stellar evolution theory. Compar- 
ison with the observations is straightforward and reflects 
directly the accuracy of the physics and the treatment of 
the (self-consistent) boundary conditions involved in the 
calculations. Any VLMS model including an adjustable 
parameter or a grey atmosphere approximation in order 
to match the observations would reflect shortcomings in 
the theory and thus would lead to unreliable results. From 
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this point of view, it is important to stress that, although 

agreement with the observation is a necessary condition 
to assess the validity of a stellar model, it is not a suffi- 
cient condition. This latter requires the assessment of the 
accuracy of the input physics and requires parameter- free, 
self- consistent calculations. 

The first generation of the present models (Baraffe et 
al. 1995) were based on the "Base" grid of model atmo- 
spheres of AUard and Hauschildt (1995; AH95). These 
models improved significantly the comparison with the 
observed Pop I and Pop II M-dwarfs sequences (Monet 
et al. 1992), down to the bottom of the main-sequence, 
w.r.t. previous models (Baraffe et al. 1995). The AH95 
models have been improved recently by including (i) a 
pressure broadening treatment of the lines, (ii) more 
complete molecular line lists, and (iii) by extending the 
Opacity- Sampling technique to the treatment of the molec- 
ular line absorption coefficients. This yields the so-called 
"NextGen" models ( Allard and Hauschildt 1997). The 
"NextGcn" synthetic spectra and colors have been com- 
pared with Population-I M-dwarf observations by Jones 
et al. (1995, 1996), Schweizer et al. (1996), Leggett et al. 
(1996), Viti et al. (1996), and have been used for the anal- 
ysis of the brown dwarf G1229B by Allard et al. (1996). 
The present work is the first application of the "NextGen" 
model atmospheres to observations and evolutionary cal- 
culations of metal-poor populations. 

A first (preliminary) set of the present improved stel- 
lar models has been shown to reproduce accurately the 
observed mass-magnitude (Henry & Mc Carthy 1993) and 
mass-spectral type (Kirkpatrick & Mc Carthy 1994) rela- 
tionships for solar-mctallicity LMS down to the bottom 
of the main sequence (Chabrier, Baraffe and Plez 1996; 
Baraffe & Chabrier 1996). The aim of the present work is 
to extend the comparison between theory and observation 
to metal-depleted abundances, characteristic of old disk 
and halo populations. 

For this purpose, we have calculated evolutionary mod- 
els for masses m < 0.8 Mq down to the hydrogen-burning 
limit, based on the "NextGen" model atmospheres, for dif- 
ferent low metal abundances characteristic of old globular 
clusters, namely [M/H] = -2, -1.5, -1.3 and -1.0. This 
low-metal grid represents a first step towards a complete 
generation of models covering solar metallicitics. A major 
source of absorption in VLMS arises from the presence 
of TiO and H2O molecular bands in the optical and the 
infrared, respectively. Although tremendous progress has 
been accomplished over the past few years in this type of 
calculations, some uncertainty remains in the absorption 
coefficients of these molecules, which affect not only the 
spectrum, i.e. the colors, but also the profile of the at- 
mosphere, and thus the evolution of late-type M dwarfs 
(see Allard et al. 1997a for a review of VLM stellar atmo- 
spheres modeling). Furthermore, the onset of grain forma- 
tion, not included in the present models, may also affect 
the spectroscopic and structural properties of metal-rich 



VLMS below %n k. 2600 K (Tsuji et al. 1996a,b; Allard 

et al. 1996). We expect these shortcomings in the models 
to be of decreasing importance with decreasing metallic- 
ity. Strong double- metal bands {TiO, VO), for example, 
are noticeably weaker in the spectra of metal-poor sub- 
dwarf (see e.g. Leggett 1992; Dahn et al. 1995), whereas 
they would dominate the spectrum of a solar metallicity 
object at this temperature. For this reason it is important 
to first examine the accuracy of the models for low-metal 
abundances and its evolution with increasing metallicity. 
The observed MS of the three aforc-mcntioncd HST GCs 
represent a unique possibility to conduct such a project. 

In some cases, comparison is made with the first gener- 
ation of models (Baraffe et al. 1995), based on the "Base" 
model atmospheres. This will illustrate the effect of the 
most recent progress accomplished in the treatment of 
LMS atmospheres. All the basic models assiimc a mixing 
length parameter, both in the interior and in the atmo- 
sphere, a = I /Hp = 1.0. Since the value of the mixing 
length parameter is likely to depend to some extent on 
the metallicity, there is no reason for the value of a in 
globular clusters to be the same as for the Sun. Although 
the choice of a is inconsequential for fully convcctive stars, 
it will start bearing consequences on the models as soon 
as a radiative core starts to develop in the interior. For 
metallicitics [M/H] < — 1, this occurs for masses m > 0.4 
Mq (see Chabrier & Baraffe 1997 for details). In order 
to examine the effect of the mixing length parameter for 
masses above this limit, we have carried out some calcii- 
lations with a = 2, both in the interior and in the atmo- 
sphere (see §3.2 below). All the models assume a solar-mix 
abundance (Grevesse & Noels 1993). As discussed in the 
previous section, oxygen-enhancement in metal-poor stars 
is taken into account by using the eqn.(l) as the correspon- 
dence law between solar-mix and a-enhanccd abundances. 
In order to test the accuracy of this procedure, a limited 
set of calculations with the exact a-enriched mixture have 
also been performed (§3.4). 

The adopted helium fraction is y = 0.25 for all metal- 
licitics. The effect of helium abundance will also be exam- 
ined below. 

All models have been followed from the initial deu- 
terium burning phase to a maximum age of 15 10^ yr. 

Before comparing the results with observations, we 
first examine some intrinsic properties of the models. 
Comparison is also made with other recent LMS models 
for metal-poor stars, namely D'Antona & Mazzitclli (1996; 
DM96) and the Teramo group of Dr. Castellani (Alexan- 
der et al. 1997). Although the EOS used by the Teramo 
group is the same as in the present calculations (SCVH), 
both groups use grey model atmospheres based on a T(r) 
relationship, and match the interior and the atmosphere 
structures cither at an optical depth r = 2/3 (DM96) or 
at the onset of convection (Teramo group). 
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3.1. Mass- effective temperature relationship 

Figure 1 shows the mass-effective temperature relation- 
ship for two mctallicities, [M/H] = -1.5 (Z = 6 x lO^'^) 
(soUd hne) and [M/H] = -1.0 = 2 x IQ-^) (dash- 
dot line). We note the two well established changes in 
the slope at Tcfr ~ 4500 K and 3800 K, respectively, for 
these metallicities. The first one corresponds to the onset 
of molecular formation in the atmosphere, and the related 
changes in the opacity (see e.g. Copcland, Jensen & Jor- 
gensen 1970; Kroupa, Tout & Gilmore 1990). The second 
reflects the overwhelming importance of electron degener- 
acy in the stellar interior near the hydrogen burning limit 
(DM96; Chabrier & Baraffe 1997). As expected, the en- 
hanced pressures of lower metallicity mixtures yield larger 
effective temperatures for a given mass (see Chabrier & 
Baraffe, 1997). 
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Fig. 1. Mass - effective temperature relation for [M / H]—-1.5 
(solid lino) and [M/H]^-1.0 (dash-dot line). The full circles 
denote our calculations done with the Krisha-Swamy grey-like 
treatment (see text) for [M/H]=-1.5 while the dashed line show 
the results of the Teramo group (Alexander et al., 1997) for the 
same metallicity. The dotted line shows the results of D'Antona 
& MazzitelU (1996; DM96) for [M/H]=-1.0. 

Also shown on the figure are comparisons with DM96 
and the Teramo group, for comparable metallicities. The 
DM96 models (dotted line) overestimate substantially the 
effective temperature, by up to ^ 200 K below 0.5 Mq. 
This will yield substantially blue-shifted sequences in the 
CMDs. This is a direct consequence of using a T(t) pre- 
scription (and a different EOS, and adiabatic gradient). 



which implies a grey atmosphere and radiative equilib- 
rium, two conditions which become invalid under low- 
mass star characteristic conditions, as demonstrated in 
previous work (Chabrier, Baraffe & Plez 1996; Chabrier 
& Baraffe 1997). The Krishna-Swamy (1966,KS) formula 
used by the Teramo group is known to slightly improve 
this shortcoming. As shown in the figure, we reproduce 
the Teramo results (dashed line) when using the KS T(t) 
relationship (filled circles). Though smaller, the depar- 
ture of the mass-Teff relationship in that case is still 
significant, in particular near the brown dwarf domain, 
where it becomes too steep and predicts too large hydro- 
gen burning minimum masses (see Chabrier & Baraffe, 
1997). A detailed examination of these different prescrip- 
tions is given in Chabrier & Baraffe (1997). These com- 
parisons show that, even though some improvement can 
be reached w.r.t. the basic Eddington approximation, any 
atmosphere-interior boundary treatment based on a pre- 
scribed T(t) relation yields incorrect, substantially over- 
estimated, effective temperatures for Tcs ^ 4500 K, i.e. 
TO ^ 0.5 Mq. Conversely, they underestimate the mass 
for a given temperature (luminosity) below this limit, i.e. 
at the bottom of the MS. This illustrates the unreliabil- 
ity of any grey-like treatment to describe accurately the 
non-grey effects and the presence of convection in opti- 
cally thin layers of VLM stars. This bears important con- 
sequences for the cooling history and thus the evolution 
in general, and for the mass-luminosity relationship and 
thus the mass calibration in particular. 

3.2. Effect of the mixing length 

In this section, we examine the effect of the mixing length 
parameter a = I /Hp on the stellar models, for a given 
metallicity. We must distinguish the effect of this param- 
eter in the stellar interior and in the atmosphere. Both 
effects reflect the importance of convective transport in 
two completely different regions of the star, namely the 
envelope and the photosphere, and thus do not necessar- 
ily bear the same consequences on the models. 

We first examine the effect of the mixing length param- 
eter in the atmosphere, by conducting evolutionary calcu- 
lations with non-grey model atmospheres for [M/H] =-1.5, 
calculated with a = 1 and 2, while the value is kept un- 
changed in the interior. We have selected a range of effec- 
tive temperatures and gravities (T^ff = 4000 - 5800 K, log 
g = 4.5 - 5) corresponding to masses m = 0.3 — 0.7 Mq. 
The atmosphere profiles corresponding to both situations 
are shown in Fig. 2, where the onset of convection and 
the location of the optical depth r = 1 are indicated. As 
seen on the figure, the atmosphere profile is rather insen- 
sitive to a variation of a in the optically-thin region, as 
expected from the rather shallow convection zone in the 
atmosphere at this metallicity and effective temperature. 
The main consequence is that the spectrum, and thus the 
colors and magnitudes at a given T^ff, remain almost un- 
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affected (less than 0.04 mag in the afore- mentioned range). 

Below 3000 K, the effect of the mixing length in the 
atmosphere is inconsequential (see also Brett f995). The 
evolutionary models calculated with both sets of model at- 
mospheres in the selected range of effective temperatures 
differ by less than 50 K in T^g. The effect of the mixing 
length in the atmosphere thus bears no consequence on 
the evolution. 
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Fig. 2. Atmosphere profiles for [M/i/]=-1.5 and a mixing 

length I — Hp (solid line) and I — 2 Hp (dash lino), respec- 
tively, for selected effective temperatures. The open circles in- 
dicate the location of t = 1 and the crosses show the onset of 
convection. 

We now examine the effect of convection in the stellar 
interior, by conducting calculations with a = 1 and 2 in 
the interior, keeping the fiducial value a = 1 in the atmo- 
sphere. This is illustrated in Figure 3a by comparing the 
solid line (and empty circles), which corresponds to a = 1 
in the interior, with the dashed line (triangles), which cor- 
responds to a = 2, for t = 10 Gyr. As expected the effect 
is null or negligible for convection-dominated interiors, i.e. 
m ^ 0.6 A/q. Above this limit, a larger mixing parame- 
ter, i.e. a more efficient convection, yields slightly hotter 
(bluer) models (D'Antona & Mazzitelli 1994; Chabrier & 
Baraffe, 1997), about 200-300 K for Toff = 5800 - 6500 
K, i.e. m = 0.7 — 0.8 Mq, for the metallicity of interest, 
[M/H] = -1.5. 

In summary, the effect of the variation of the mixing 
length parameter on the stellar models remains weak for 
the lower MS and is essentially affected by the value in the 
stellar interior. A variation of a = 1 to 2 yields variations 



ATeff/Teff 2 — 5% on the effective temperature in the 
upper MS, for masses above ^ 0.6 Mq, while the lumi- 
nosity remains essentially unaffected. We will examine in 
section 4.1 whether comparison with the observed MS of 
GCs can help calibrating this value. 

3.S. Effect of the age 

The effect of the age in a color-magnitude diagram is 
shown on Figure 3a, for [M/H] = — 1.5 for ages t = 10 
and 15 Gyr. The time required to reach the zero-age main 
sequence is much smaller than the age of GCs over the 
entire stellar mass range so that all hydrogen burning ob- 
jects do lie on the MS and the bottom of the GC MS 
is unaffected by age variation. On the other hand, metal- 
poor stars are significantly hotter and more luminous than 
their more metal-rich counterparts. Therefore, for a given 
mass, they burn more rapidly hydrogen in their core, and 
thus evolve more quickly off the MS. The more massive 
(i.e. hotter) the star, the larger the effect. This defines 
the turn-off point, i.e. the top of the MS. For fixed metal- 
licity, the turn-off point will thus be reached for lower 
masses, i.e. fainter absolute magnitudes, as the age in- 
creases. This is illustrated in Figure 3a where we compare 
the 10 Gyr and 15 Gyr isochrones calculated for different 
mixing length parameters, for the same metallicity. While 
the position in the HR diagram of masses below ~ 0.6 Mq 
remains unchanged, larger masses become substantially 
bluer and more luminous with age as they transform more 
central hydrogen into helium. The resulting increase of 
molecular weight yields further contraction and heating 
of the central layers, and thus an increase of the radia- 
tive flux {Frad oc T^) and the luminosity. Eventually they 
evolve off the MS for a (turn-off )-mass mro ~ 0.8 Mq for 
[M/H] = -1.5 and -2 at t=10 Gyr and mro 0.75 Mq 
at 15 Gyr. As shown on the figure, the upper MS is sensi- 
tive to age and mixing length variations. The calibration of 
the mixing length parameter on the turn-off point is then 
altered by age uncertainties. The precise determination of 
the turn-off point thus requires more detailed calculations, 
which are out the scope of the present study. 

3.4- Effect of helium and a-elements 

As mentioned above, the fiducial calculations have been 
conducted with a helium abundance Y = 0.25. Figure 3b 
compares these results (solid line; circles), for [M/H] = 
— 1.5, with calculations done with Y = 0.23 (-I-). The cor- 
responding variations of the temperature (color) and the 
luminosity (magnitude) are negligible over the entire mass 
range, except near the turn-off point. 

In order to examine the accuracy of the Ryan & Nor- 
ris prescription (eqn.(l)), we have conducted calculations 
using model atmospheres computed with an a-element 
abundance enrichment of [a/Fe] = 0.3 for a [Fe/H] = 
-1.3 {[a/H] = -1.0) mixture. Since the spectra of VLMS 
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Fig. 3. (a) Effect of the age and of the mixing length in 
the stellar interior for [M/H]=-1.5. The models display the 
t = 10 Gyr isochrone when using a mixing length / = Hp 
(solid line; empty circles) and I = 2 Hp (dashed line; trian- 
gles). The dotted line (+) shows the 15 Gyr isochrone calcu- 
lated with / = Hp, whereas the dash-dotted curve {X) corre- 
sponds to t=15 Gyrs and / = 2 Hp. The masses indicated on 
the figure correspond to the 10 Gyrs isochrones (open circles 
and triangles). The upper mass corresponding to the 15 Gyrs 
isochrones is 0.75 Mq. (b) Effect of the helium fraction Y, for 
[M/H] = —1.5, with Y = 0.25 (solid curve, open circles) and 
Y = 0.23 {+). 



near the bottom of the MS are governed by a-elements, we 

expect this mixture to yield equivalent results to those ob- 
tained with the scaled-solar abundance mixture of [M/H] = 
-1.0. This test will also determine whether neglecting a rel- 
ative under-abundancc of Fe and other non-a elements in 
the scaled-solar abundances does affect or not the shape 
of the MS. We have computed model atmospheres with 
the afore-mentioned modified abundances for several tem- 
peratures and gravities in the range Teff = 3000 — 5000 
K and logg = 4.5 — 5.0. We find that the atmosphere 
density and temperature profiles and the colors are es- 
sentially undistinguishable from the ones calculated with 
[a/Fe] = 0.0, [M/H]{= [Fe/H]) = -1.0, i.e. a scaled 
solar-mix, throughout the entire atmosphere. 

We have also examined the effect of a-enriched abun- 
dances in the mienor opacities (OPAL), thanks to appro- 
priate opacities kindly provided by F. Rogers for [M/H]=- 
1.0 and [a/Fe] = 0.3. Here too we find that the effect is 
negligible and does not modify noticeably the isochrones. 
This assesses the validity of the Ryan & Norris procedure 
to take into account the a-element enrichment for metal- 
poor stars, and demonstrates the consistency of our pre- 
scription, i.e. comparing solar- mix models with [a/H] = 
[M/H] = -1.0 to observations with [Fe/H] = -1.3. Con- 
versely, this demonstrates the inconsistency of compar- 
ing observations and solar-mix models for these stars with 
[Fe/H] — [M/H], as done sometimes in the literature. 

The theoretical characteristics of the present mod- 
els, effective temperature, luminosity, gravity, bolomet- 
ric magnitude and magnitudes in VRIJHK for several 
metallicities and an age t = 10 Gyrs are given in Tables 
II- V. 

4. Results cind discussion 

4-1- Comparison with globular cluster main sequences 

Figures 4-6 show the main sequence CMDs of the three 
GCs mentioned in §2. As mentioned previously, the CMDs 
are shown in the M555 (~ V), M814 (~ /), and Meoe 
filters in the WFPC2 Flight system to avoid possible errors 
due to uncertain photometric conversion into the standard 
Johnson-Cousins system. In all cases, the main sequence is 
well defined down to J « 25, V ~ 27. Below this limit, the 
observed MS dissolves into the field and it becomes quite 
difficult to distinguish the cluster-MS from field stars. 

The theoretical MS for the appropriate metallicity, as 
described in §2, are superimposed to the observations in 
each figure. For comparison in the Flight system, observed 
magnitudes have been dereddened with the corrections de- 
rived from the synthetic spectra (cf. Table I). The first 
striking result is the excellent agreement between theory 
and observation for the three clusters, spanning a fairly 
large metallicity range from strongly metal-depleted abun- 
dances (M15) to a tenth of solar metallicity (ojCen). In 
particular, the changes of the slope in the observed MS are 
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perfectly reproduced by the models, for the correct magni- 
tude, color and metallicity. Since these changes stem from 
the very physical properties of the stellar interior and at- 
mosphere, as discussed in the previous section, the present 
qualitative and quantitative agreement assesses the accu- 
racy of the physical inputs in the present theory. The 
masses corresponding to these changes are indicated on 
the curves and their effective temperatures are given in 
Tables II - V, for the various metallicities. 

Let us now consider each cluster in turn in order of 
increasing metallicity. 

• M15 (Fig. 4a, b) : This is the most metal-depleted 
HST cluster presently observed, with [Fe/H] ~ —2.26 
(Djorgovski 1993), i.e. [Al/H] ^ —2.0 within the prescrip- 
tion adopted in the present paper. The MS of M15, as ob- 
served by DeMarchi & Paresce (1995), is shown in Figure 
4a, with the distance modulus adopted by these authors. 
Comparison is made with our models with [M/H] = —1.5 
(solid line) and -2.0 (dashed line). Although both mod- 
els merge in the upper MS, metallicity effect is clearly re- 
flected in the lower MS where the peak of the spectral dis- 
tribution falls in the V bandpass, i.e. below m < 0.6 M©, 
Mgoe - M814 > 0.6, and Toff < 5000 K. As seen on the 
figure, the agreement with the observations over the whole 
MS is excellent for [M/H] = —2 (the observationally- 
determined metallicity) although the models appear to 
be slightly too red by ~ 0.05 mag or over luminous by 
~ 0.2 mag, in the upper MS. Even though this discrepancy 
is within the observational error bars in mgoe and m8i4, 
which range from ±0.03 mag for the upper MS to ±0.17 
mag for the lower MS (DeMarchi, private communication), 
we will examine the possibilities for this disagreement to 
be real. A slightly lower metallicity would leave the upper 
MS almost unchanged, as seen from the two theoretical 
sequences displayed on the figure. Although the offset in 
the upper MS might be compensated by a slightly larger 
mixing length (see §3.2), the fact that it is quite constant 
along the entire sequence — whereas the mixing length 
only affects the upper MS (see §3.2) — rather suggests 
either an overestimated reddening E{B — V) (by ^ 0.05 
mag) or an underestimated distance modulus. Indeed, as- 
suming (m — M)si4 = 15.5, i.e. ±0.2 mag w.r.t. the value 
quoted by DeMarchi & Paresce, or a reddening fainter by 
~ 0.05 mag would bring theory and observation into per- 
fect agreement. A 0.2 mag error on the distance modulus 
corresponds to 1 kpc, i.e. ^ 10% of the canonical value 
used presently (Djorgovski 1993). Such an error cannot be 
excluded for this remote cluster {D w 10.5 kpc). 

Note that we exclude an artificial offset in the cali- 
bration of the model magnitudes as source of discrepancy, 
since we use the magnitude zero-points kindly provided 
by G. DeMarchi. However, a calibration problem of the 
data in the F606W" filter cannot be excluded, since a dis- 
agreement appears as well in the comparison of NGC6397 
observed by the same group in the same filters. The dis- 
agreement however vanishes for the same cluster in M505 



and M814 and for wCen observed by Elson et al. (1995) in 

Meoe and M814 (see below). The latter agreement seems 
to reject an intrinsic problem of the atmosphere models in 
the spectral region covered by the Mgoe passband. 

The data of DeMarchi & Paresce have been trans- 
formed into the standard V—I Johnson- Cousins system by 
Santiago et al. (1996) who used the dereddening correction 
quoted by the former authors (see Table I) and the syn- 
thetic transformations of Holtzman et al. (1995). Figure 
4b compares the M15 MS fit derived by Santiago et al. to 
the present models in the same Johnson-Cousins color sys- 
tem, where the filter transmissions of Bessell (1990) have 
been applied to the models. The agreement is similar to 
that obtained in the previous Flight system, well within 
the error bars of the photometric meoe-to-Vj transforma- 
tion, with the same slight offset in color. This comparison 
assesses the validity of the photometric synthetic transfor- 
mations of Holtzman et al. (1995) in the present filters. 

As shown in the figure and in Table II, the faintest 
observed stars on the MS, M814 10 or My - 11.5, 
correspond to a mass m ~ 0.15 M©, still well above the 
hydrogen burning limit. This latter, m = 0.083 for 
[M/H] = —2, is expected to correspond to M814 ^ 14.6, 
Meoe - M814 ~ 4.3 or My ~ 21 and (F - /) ~ 6.5 (Table 
II). At the bright end, the limit of the observations corre- 
sponds to the turn-off point M814 ~ Mj ^ 3 (DeMarchi 
& Paresce 1995) which corresponds to ~ 0.8 ± 0.05 Mq 
depending on the age. Note that the Johnson-Cousins fit 
given by Santiago et al. (1996) starts at My ^ 6, which 
corresponds to a mass m ~ 0.6 — 0.7Mq and does not 
include the upper MS. 

• iVG'C6397 (Fig. 5a,b,c) : This cluster has been ob- 
served separately by Paresce et al. (1995) and Cool et al. 
(1996) in different photometric filters, thus allowing com- 
parison in the three WFPC2 filters mentioned previously. 
Figure 5a shows the comparison of the observations and 
the models in M814 vs Meoe - M814 for [M/H] = -1.5 
and —2. The best agreement is obtained for the observed 
metallicity [M/H] = —1.5, although admitedly the the- 
oretical sequence lies near the blue edge of the low MS 
(M814 <; 8). We note however that the error bars for 
this cluster range from ±0.02 for the brightest stars to 
± 0.20 mag for the lower MS (cf. Paresce et al. 1995). The 
agreement displayed in Fig. 5a is therefore well within the 
error bars. Much better agreement is found in M555 vs 
— Msi4, as shown on Figure 5b, although an off- 
set now appears for the intermediate part of the sequence 
(-M555 ~ 8), ~ 0.05 mag in color and ~ 0.3 mag in mag- 
nitude. As shown on Figure 3a, this part of the MS is 
insensitive to the mixing length, and the disagreement 
would thus not be solved by a larger mixing length. In 
the same vein, models with a-enriched abundances would 
yield similar results, as discussed in §3.4. Models with a 
substantially lower metallicity would fail reproducing the 
shape of the MS both in M555 (cf. Fig. 5b) and M814 (see 
Fig. 5a). 
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Fig. 4. (a) CMD for M15. The data in the 606 and 814 filters are from DeMarchi & Paresce (1995). The data are dereddened 

with the corrections calculated in the present study, as given in Table I. Solid line (o): [M /H]—-1.5; dashed line (+): [M / H]—-2.0. 
(b) Same as in Fig. 4a in the Johnson- Cousins system with the fit of Santiago et al. (1996) for the cluster (thick dash-dot). 



Both figures 5a and 5b are displayed with the distance 

moduhis quoted by Cool ct al. (1996), i.e. (m — M)q = 
11.7. The choice of the value quoted by Paresce et al. 
(1995 ), (m - M)o = 11.9, brings the tracks for [M/H] = 
— 1.5 exactly on the observed sequence for the M555 vs 
-^555 — Msi4 CMD, as illustrated on Figure 5c. Such an 
undetermination in the distance modulus of NGC6397 
corresponds to a difference of 200 pc, i.e. about 10%. In 
the Mgoe and Msi4, filters, the models are slightly too 
blue, with an offset in color ;^0.05 mag, which can well 
stem from the same calibration problem in F606W men- 
tioned previously for the cluster M15. The predicted se- 



quence with the metallicity [M/H]=-1.5 remains however 
well within the error bars of the observed MS. 

It is interesting to analyse the agreement obtained by 
the Teramo group models for this cluster. We first note 
that their models, while based on a solar-mix, correspond 
to [M/H] = [Fe/H] rather than [M/H] = [O/H]. As dis- 
cussed in the previous section, this leads to inconsistent 
comparisons, and is reflected by the fact that: i) they fit 
the MS with a rather low value of the metallicity for this 
cluster (-1.9 < [M/H] < -1.6), and u) they adopt a 
reddening correction E{V — I) = 0.19 which differs sig- 
nificantly from the value prescribed by Cool et al. (1996) 
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Fig. 5. (a) CMD for NGC6397 in the 606 and 814 filters. The data are from DeMarchi & Paresce (1995). The data are 
dereddened with the corrections calculated in Tabic I. Solid line (o): [M / H]=-1.5; dashed line (+): [M / H]=-2.0. The distance 
modulus is the one quoted by Cool et al. (1995), i.e. (m - M)o=11.7. (b) CMD for NGC6397 in the 555 and 814 filters. The 
data are from Cool et al. (1995). The data are dereddened with the extinctions calculated in Table I. Dot-dashed line (•): stellar 
models based on the "Base" model atmospheres (Baraffe et al. 1995) for [M/H]=-1.5; solid line (o): present stellar models based 
on the "NextGen" model atmospheres for [M/H]=-1.5; dashed line (+): [M/ H]=-2.0 with the NextGcn model atmospheres. 
The distance modulus is the same as in Figure 5a, i.e. (m — M)o=11.7. (c) Same as figure 5b, with the distance modulus quoted 
by DeMaxchi & Paresce (1995), i.e. (m — M)o=11.9 and the same extinction as in Figure 5b. Solid line (o): [M/H]=-1.5; dashed 
line (-I-): [M/H]=-2.0. A comparison with the models of the Teramo group for [M/H]=-1.5 is shown (dash-dot) 



and by the present calculations (see Table I). Using the 
correct extinction would redshift significantly their models 
for [M/H]=-1.5 w.r.t. the observations, as shown on Fig. 
5c (dash-dotted line). 

Note that for this cluster, the observations are 
still above the hydrogen-burning hmit. The HBMM ^ 
0.083 Mq corresponds to My ~ 19.5, Mj ~ 13.9 for 
[M/H] = —1.5, whereas the bottom of the observed MS, 
My ~ 13, Ml ~ 11, corresponds to to = 0.11 Mq. The 
observations of Cool et al. (1996) extend to brighter mag- 
nitudes and reach the turn-off point, My ~ 4, m ~ 
0.8 ±0.05Mq. 

A limited set of the first generation of the present 
models at [M/H] = —1.5, based on the "Base" model 
atmospheres (Baraffe et al. 1995), is also shown on Fig- 
ure 5b for m < 0.2 (full circles, dash-dot line). These 
models have a slightly different trend w.r.t. the present 
ones. This stems from a general overestimation of the 
molecular blanketing, the main source of absorption in the 
coolest VLMS, because of the straight-mean approxima- 



tion in the "Base" models, and thus an underestimation 
of the flux in the V-band ( Chabrier et al. 1996; AUard & 
Hauschildt 1997). The better agreement with the present 
models clearly illustrates the recent improvements in the 
treatment of the molecular opacities, especially for TiO 
(AH97), which strongly affect the atmosphere profile, and 
thus the evolution. However, the agreement between the 
previous models and observations was already quite sat- 
isfactory (~ 0.2 mag) and for the first time reproduced 
accurately the bottom of the observed MS for metal-poor 
LMS. 

• ui — Cen (Fig. 6) : this cluster has been observed 
by Elson et al. (1995) in the F606W and F814W filters. 
The observed metallicity usually quoted for this cluster is 

[Fe/H] 1.6, i.e. [M/H] 1.3 (cf Elson et al. 1995 

and references therein), although a value [Fe/H] ~ —1.2, 
i.e. [M/H] ^ —1.0 has been suggested recently (Norris et 
al., 1996). The data presented in Figure 6 have been re- 
calibrated (Elson, private communication), since the cal- 
ibrations of Holtzmann et al. (1995) used in Elson et al. 
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Fig. 6. (a) CMD for oi-Cen. The data in 606 and 814 filters are from Elson et al. (1995). The models correspond to [M/H]=-1.3 
(solid line, open circle) and [M/H]=-1.0 (dashed line, -|-) with I = Hp. Pull circles correspond to [M/H]=-1.0 and I = 2Hp 

(b) CMD for a; — Cen in the Johnson-Cousins system. The data from Elson et al. (1995) have been recalibrated with updated 
transformations (Elson, private communication). Same models as in (a). 



(1995) were preliminary at this time. The zero points for 

both HST filters arc updated (cf. Holtzmann ct al. 1995), 
as well as the transformation into the Johnson- Cousins 
system. Here again, the match between theory and obser- 
vation is almost perfect, as shown in Figure 6, in both the 
HST instrumental and the Jonhson- Cousins systems. Note 
that both systems predict the same masses for the upper 
MS (to - O.75M0) and the lower MS (m - O.ISM©), 
well within the error bars due to the photometric con- 
version from rriQOQ into Vj (Holtzman et al. 1995). Once 
again, this assesses the accuracy of the present photomet- 
ric conversions of HST magnitudes into Johnson-Cousins 
magnitudes. 



We used the same distance modulus as Elson et al. 
(1995), i.e. (to-M)o = 13.77. Adopting the slightly higher 
value 13.92 quoted by Santiago et al. (1996) will shift the 
blue edge of the observed MS on the theoretical isochrone 
with [M/H]=-1.3, yielding a less good agreement. Our red- 
dening corrections (cf. Table 1) are in excellent agreement 
with the values quoted by Elson et al. (1995). Calculations 
with a mixing length I = 2Hp for the [M/H]=-l isochrone 
are also shown in Figure 6. Note that the observations re- 
main compatible with a value 1 ^ I /Hp ^ 2. Since the 
same remark applies to the lower-metallicity chisters anal- 
ysed previously, it would be hazardous to try to derive ro- 
bust conclusions about a possible dependence of the mix- 
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ing length on metallicity. The observed width of the lower 

MS docs not allow a precise determination of the metal- 
licity of the cluster from the theoretical isochrones. The 
shape of the MS is well reproduced with both [M/H]=- 
1.3 and -1. This result is consistent with the spread in 
metallicity determined by Norris et al. (1996). 

The properties of the models for the afore-mentioned 
different metallicities, in particular the mass-color- 
magnitude relationships, are displayed in Tables II-V. The 
lowest mass in the table corresponds to the hydrogen- 
burning minimum mass (HBMM) for each metallicity. For 
an age t = 10 Gyr and metallicities [M/H] < —1.0, the 
magnitude of the most massive brown dwarf is ^2 mag in 
/ and K, ~3 mag in V and ~ 1.5 mag in J and H fainter 
than the one corresponding to the HBMM. This sets the 
scale of the detection limit for the search for brown dwarfs 
in globular clusters with future space-based observations. 

As already mentioned and clearly seen from the fig- 
ures and from Tables II - V, the more metal-depleted 
clusters have bluer main sequence colors, a consequence 
of the increasing effective temperature at a given mass 
with decreasing metallicity, since the same optical depth 
corresponds to a denser layer in an increasingly transpar- 
ent atmosphere (see e.g. AH95; Chabrier & Baraffe 1997). 
Metallicity effects are most apparent in the intermediate 
MS, i.e. TO ^ 0.5 Mq and Tgs <, 4500 K, where the peak of 
the spectral distribution falls near the V bandpass. This 
stems essentially from the increasing TiO-opacity, which 
absorbs mainly in the optical and reddens the {V — I) 
color. Thus, for a given Tes, a metal poor object will ap- 
pear bluer than the more metal-rich coimterpart. This ef- 
fect is strengthened at fixed mass by the fact that the lower 
the metallicity, the hotter (bluer) the star. 

4-2. Color-magnitude diagram in the near infrared 

The observations of GCs in near IR colors will soon be 
possible with the next generation of HST observations, i.e. 
the NICMOS camera, and in a more remote future, with 
the european Very-Large-Telescope (VLT). The NICMOS 
filters include Wide (W), Medium (M) and Narrow (N) 
bandpasses from 1.1 to 2.4 /xm. In figures 7a-b we show 
our models at t=10 Gyrs for [M/H]=-2, -1.5 and -1.0 in 
the NICMOS wide filters FllOW, F160W and F187W. 
For comparison, the same isochrones are displayed in the 
J (~ HOW) and H (- 160W) magnitudes (cf. Fig. 7a, 
dotted curves), defined in the CIT system (Leggett 1992). 
The masses listed in Table II-V are indicated by the signs 
(-1- and circles) on the curves (except 0.13 Mq excluded 
for sake of clarity). We note the ongoing competition, in 
these infrared colors, between the reddening due to the 
decreasing temperature and increasing metallic molecu- 
lar absorption in the optical and the increasing collision- 
induced absorption of H2 in the infrared (cf. Saumon et 
al., 1994; AUard & Hauschildt, 1995) which shifts back the 



flux to shorter wavelengths. This leads to quasi-constant 

color sequences from ~ 0.5 to ^ 0.1 Mq, corresponding to 
Teff ~ 4500 K, for which H2 becomes stable in the atmo- 
sphere, to ~ 3500 K, as predicted also for zero-metallicity 
(Saumon et al., 1994). Below this limit, molecular hy- 
drogen becomes dominant, the density keeps increasing 
( cf. Chabrier & Baraffe, 1997) and H2 CIA-absorption 
becomes the dominant effect {kcia ~ KH2-H2 P^h^ ^"^^^ 
order, see Guillot et al., 1994). This causes the blue loop at 
the very bottom of the MS in IR colors, as seen in Figures 
7, whereas the optical colors redden monotonically with 
decreasing mass, as shown in the tables. The blue loop be- 
comes more dramatic with decreasing metallicities, since 
the lower the metallicity the denser the atmosphere. This 
general trend is very similar in all NICMOS filters cover- 
ing the above-mentioned wavelength range and is clearly 
a photometric signature of the stellar to sub-stellar tran- 
sition, whose physical source is the large increase of the 
density in this region and the ongoing H2 molecular re- 
combination and collision-induced absorption. The limit 
magnitudes required to reach the very bottom of the MS 
are ~ 14 for [M/H]=-2, ~ 12.5 for [M/H]=-1.5 and ~ 12 
for [M/H]=-l, and are essentially the same for all NIC- 
MOS filters from FllO to F240. The more massive brown 
dwarfs will be about 2 magnitudes fainter. 



4-3. Comparison with the halo suhdwarf sequence 

Figure 8 displays different observations of halo field stars 
in the standard Johnson-Cousins system. The filled cir- 
cles are the sM&rf7j;arf sequence of Monet et al. (1992), the 
crosses the more complete subdwarf sequence of Dahn et 
al. (1995), and the triangles correspond to a sub-sample 
of Lcggett's (1992) halo stars. The halo classification was 
determined photometrically (Leggett 1992) and kinemati- 
cally : the stars in the three afore- mentioned samples have 
tangential velocities \Vtan\ > 220 km.s"-'^ (Monet et al.), 
>160 km.s~^ (Dahn et al.) and >180 km.s~^ (present 
sub-sample of Leggett). All these observations appear to 
be fairly consistent, the Monet et al. sample representing 
the most extreme halo fraction of the Dahn et al. sample, 
the Leggett's sample containing only a few genuine halo 
stars. We stress that the linear fit proposed by Leggett 
(1992) is a rather poor representation of the distribution 
of the true halo objects and is strongly misleading. We 
also emphasize that linear fits are not correct to fit LMS 
sequences in the HR diagram, since they do not repro- 
duce the wavy behaviour of the sequences, wich reflects 
intrinsic physical properties of these stars, as discussed in 
§3.1. This non linearity has already been stressed for the 
mass-luminosity (D'Antona & Mazzitelli, 1994; Chabrier 
et al., 1996) and the mass-spectral class relations (Baraffe 
& Chabrier, 1996). We first note the important spread in 
color, over 1 mag, which reflects the large spread in metal- 
licity in the sample. 
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Fig. 7. (a) Isochroncs at t=10 Gyrs in the NICMOS filters (solid curves) FllOW and F160W for three different mctallicity 
[M/H]=-2, -1.5 and -1 (from left to right). The dotted curves correspond to the same isochrones in the CIT system Mj-(J-H). 
The signs on the curves correspond to the masses (except 0.13 Mq) tabulated in Tables II to V. (b) Same as in Fig. 7a in the 
NICMOS filters FllOW and F187W. 



The solid lines indicate the present LMS sequences for 
[M/H] = -2.0, -1.5, -1.3 and -1.0 from left to right. The 
Monet et al. sample is consistent with an average mctallic- 
ity [M/H] ~ —1.3 to —1.5. It is more hazardous to try to 
infer an average value for Leggett's sample, given the lim- 
ited number of objects and the large dispersion. The Dahn 
at al. sample clearly includes the two previous ones, with 

objects ranging from [M/H] 0.5 to [M/H] = -2.0. 

This reflects the difficulty to determine the precise origin 
of an object from its kinematic properties only. 

As seen on the figure, these field stars are fairly con- 
sistent with the theoretical sequences determined for the 
clusters, for similar metallicities. This is particularly ob- 
vious for the sequences of iVGC6397 ([M/H]=-1.5) and 
ojCen ([M/H]=-1.3), which match perfectly the Monet 
et al. average sequence. This strongly suggests, contrar- 



ily to what has been suggested by Santiago et al. (1996), 
that there is no significant difference in the structure and 
the evolution of globular star clusters and field halo stars. 
The discrepancy between the CMD of Leggett (1992) and 
those from the HST led Santiago et al. (1996) to invoke the 
possibility of a calibration problem of the HST data. How- 
ever, the agreement that we find between our theoretical 
models and the globular cluster CMDs in both the HST 
data and the Johnson-Cousins system (e.g for M15 and 
wCen), excludes this hypothesis. It rather stems form the 
large metallicity dispersion in Leggett's sample and from 
the misleading fit of this sample. 
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Fig. 8. Subdwarf halo field stars from the data of Monet et al. (1992) (full circles), Dalin et al. (1995) (+) and Leggctt 
(1992) (triangles), in the standard Jonhson-Cousins system. All the stars shown here from the Leggett's sample have tangential 
velocities \Vtan\ > 180 km.s~^. Solid lines: present models for [M/H]=-2,-1.5, -1.3 and -1 (from left to right) for t=10 Gyr. The 
empty circles correspond to the masses indicated on the figure. 



5. Conclusion 

Wc have presented extensive calculations of VLMS evo- 
lution in the range —2.0 < [M/H] < —1.0, characteristic 
metallicities for old globular star clusters and halo field 
stars. These low-metallicities minimize possible shortcom- 
ings of the model atmospheres pertaining to incomplete or 
inaccurate metallic molecular line lists and grain forma- 
tion, and provide a stepping stone towards the derivation 
of more accurate low-mass star models for solar metal- 
licity. The models are examined against available deep 
photometry color-magnitude diagrams obtained with the 
Hubble Space Telescope for three globular clusters. The 
HST CMDs for the clusters span a large range in metal- 
licity, thus providing very stringent tests for the models. 
Since the parameters characteristic of these clusters, ex- 
tinction, distance modulus and metallicity are fairly well 
defined, there is no free parameter left to bring models 
into agreement with observations. Therefore comparison 
between theory and observation reflects directly the ac- 
curacy of the theory. We stress the importance of eqn.(l) 



when comparing GC observed and theoretical CMDs. A 
first generation of the present models (Baraffe et al. 1995) 
has been used incorrectly by comparing observations at 
[Fe/H] with a model at the same value of [M/H]. 

The main conclusions of these calculations can be sum- 
marized as follows : 

• We first note the overall remarkable agreement be- 
tween the present models and the observations, within less 
than 0.1 mag, over the whole metallicity range and the en- 
tire main sequence from the turn-off to the bottom. The 
characteristic changes in the slopes of the cluster MS's 
are reproduced accurately, and assess the validity of the 
physics involved in the models. This yields an accurate cal- 
ibration of the observations, i.e. reliable mass-magnitude- 
effective temperature-age relationships. 

We also provide reddening corrections based on accu- 
rate LMS synthetic spectra. 

• Variations of the mixing length in the stellar inte- 
rior affect essentially the upper main sequence near the 
turn-off, i.e. only stars massive enough to develop a large 
radiative core (m ^ 0.7 M©). Variations of the mixing 
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length in the atmosphere is found to be inconsequential 
on cvohitionary models. There is no clear hint for a de- 
pendence of the mixing length on the metallicity and our 
results remain in agreement with all observed sequences, 
within the error bars, for 1 < Imix/Hp < 2. 

• The Ryan & Norris (1991) prescription to convert 
solar-mix abundances into oxygen-enriched mixtures char- 
acteristic of old stellar populations is accurate. Using the 
correct a-enhanced mixture (both in the atmosphere and 
in the interior opacities) leads to results identical to those 
obtained with the aforc-mcntioned scaling. Not taking 
this enrichment into account leads to inconsistent com- 
parisons. 

• We derive theoretical sequences in the filters of the 
NICMOS camera, down to the hydrogen burning limit. 
This will allow a straightforward analysis of the future 
HST observations, and provides a stringent test for the 
accuracy of the present models near the brown dwarf limit. 
This corresponds to Mj « Mh » Mk « 13 - 14 for the 
lowest metallicity examined presently, i.e. [M/H] = —1.0. 

We also predict a photometric signature of the tran- 
sition from stellar to substellar objects in the infrared, 
in terms of a severe blue loop near the very bottom of 
the MS, whereas optical colors keep reddening almost lin- 
early. This photometric signature reflects the overwhelm- 
ing absorption of molecular hydrogen in the infrared due 
to many-body collisions, and stems from the increase of 
the molecular hydrogen fraction and of the density (con- 
traction) near the stellar to substellar transition. 

• The models allow a good determination of the metal- 
licty of the observed halo field stars. A striking result is the 
large metallicity dispersion of these objects, from [M /H] 
= -2 to near solar, although they all have halo kinematic 
properties. The most extreme halo stars, represented by 
the Monet et al. (1992) sample, have a metallicity rang- 
ing from [M/H] ~ —1.0 to -2.0, with an average value 
[M/H] ~ —1.3. We find no evidence for differences in 
the sequences of halo field subdwarfs and globular star 
clusters. Both are reproduced with the same isochrones 
for similar mean metallicity, [M/H] ~ —1.3 to —1.5, i.e. 
[Fe/H] 1.6 to -1.8. 

We can now affirm that the theory of low-mass stars, 
at least for metal-depleted abundances, has reached a very 
good level of accuracy and can be used with confidence to 
analyse the observations and make reliable predictions. 
While we cannot honestly exclude the possibility of a dis- 
crepancy in the Mqqq — Msu color predicted by the models 
(by ~ 0.05 mag), it remains within the error bars due to 
the observations themselves and to the undeterminations 
of either the extinction or the distance modulus. It may 
also reflect present uncertainties in the calibration of the 
F606 filter. 

While it is unfortunate that the hydrogen burning limit 
remains unreached by the present HST observations, for 
it is masked by foreground field stars, the next cycles of 
HST observations will be able to resolve it in the near fu- 



ture. Indeed, known proper motion of the GCs will cause a 

substantial displacement of cluster stars between 1994 and 
1997 which should allow a separation of cluster from field 
stars (King 1995). Also, as mentioned above, observations 
in the IR should lead to the observation of the very bot- 
tom of - and possibly below - the MS. The present models 
provide the limit magnitudes to be reached to enter the 
brown dwarf regime. 

The present calculations represent an important im- 
provement in the description of the mechanical and ther- 
mal properties of low-mass stars, and of their photometric 
signature. This provides solid grounds to extend these cal- 
culations into the more complicated domain of solar-like 
metallicities, as will be examined in a forthcoming paper 
(AUard et al., 1997b). The assessed accuracy of the present 
models provide reliable mass-luminosity relationships for 
metal-poor stellar populations in general and for globu- 
lar clusters in particular. This allows for the first time 
the derivation of reliable mass-functions for these objects 
down to the brown dwarf limit (Chabrier and Mera, 1997). 

Tables H-V are available by anonymous ftp: 
ftp ftp.ens-lyon.fr 
user name: anonymous 
ftp > cd /pub/uscrs/CRAL/ibaraffe 
ftp > get BCAH97.models 
ftp > quit 
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Table I. Characteristics of the globular clusters. [Fe/H] is 
the observed metallicity, [M/H] is the metallicity used in the 
models, which takes into account the a-enrichment following 
the prescription of Ryan & Norris (see §2). The fourth and 
fifth columns denote the bare distance modulus and redden- 
ing repectively, as quoted by the refered observers. The last 
three columns give the extinction in the three respective fil- 
ters F555 ~ Vj, F606 and _F814 ~ /c, where J and C denote 
the standard Johnson- Cousins system. The upper rows give 
the value chosen by the refered observers whereas the lower 
rows denote the value used in the present calculations, deduced 
from the AUard & Hauschildt (1997) synthetic spectra and the 
Cardelli et al. (1989) extinction law. In all cases the interstellar 
extinction corresponds to Rv = 3.12. 



cluster 


[Fe/H] 


[M/H] 


(m-M)o 


E(B-V) 


^555 


^606 


^814 


M15^'^ 


-2.3 


-2.0 


15.10 


0.11 


0.34 


0.274 


0.144 








15.10 


0.11 


0.33 


0.30 


0.20 


NGC63973 


-1.9 


-1.5 


11.70 


0.18 


0.57 




0.34 


NGC6397'' 


-1.9 


-1.5 


11.90 


0.18 






0.30 








11.7 (11.9) 


0.18 


0.57 


0.51 


0.35 


u> Cen^ 


-1.6 


-1.3/-1.0 


13.77 


0.116 


0.36 


0.32 


0.22 








13.77 


0.116 


0.36 


0.33 


0.22 



^ Piotto et al. 1995 

^ De Marchi & Parcscc 1995 

^ Cool et al., 1996 

* Paresce et al. 1995 

^ Elson et al. 1995 



This article was processed by the author using Springer- Verlag 
M^X A&A style file L-AA version 3. 
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Table II. Physical properties and absolute magnitudes of low-mass stars for [M/H] = —2.0 and t = 10 Gyrs. The lowest mass 
corresponds to the hydrogen-burning limit. The mass m is in Mq, Teff in K and the luminosity L in Lq. The VRI magnitudes 
are in the Johnson- Cousins system (Bessell 1990) and the JHK magnitudes in the CIT system (Leggett 1992). Note that the 
bolometric magnitude corresponds to Mhoi{Q) = 4.64. 



m 


Teff 


log L 


log g 


Mbol 


Mv 


Mr 


Mi 


Mj 


Mh 


Mk 


0.083 


1779. 


-4.274 


5.58 


15.33 


20.92 


16.92 


14.41 


13.47 


13.95 


14.58 


0.085 


2342. 


-3.741 


5.53 


13.99 


17.67 


15.17 


13.11 


11.98 


12.30 


12.61 


0.090 


3005. 


-3.202 


5.45 


12.65 


14.70 


13.40 


11.96 


10.77 


10.60 


10.68 


0.100 


3395. 


-2.860 


5.37 


11.79 


13.25 


12.22 


11.17 


10.06 


9.66 


9.58 


0.110 


3565. 


-2.686 


5.32 


11.36 


12.63 


11.66 


10.76 


9.69 


9.24 


9.12 


0.130 


3749. 


-2.463 


5.26 


10.80 


11.91 


10.99 


10.21 


9.20 


8.72 


8.58 


0.150 


3852. 


-2.306 


5.21 


10.40 


11.43 


10.55 


9.81 


8.84 


8.36 


8.20 


0.200 


4004. 


-2.020 


5.12 


9.69 


10.59 


9.78 


9.09 


8.17 


7.69 


7.52 


0.300 


4172. 


-1.665 


5.01 


8.80 


9.58 


8.83 


8.19 


7.33 


6.85 


6.69 


0.400 


4304. 


-1.421 


4.94 


8.19 


8.87 


8.18 


7.58 


6.75 


6.29 


6.13 


0.500 


4623. 


-1.083 


4.83 


7.35 


7.83 


7.25 


6.72 


5.99 


5.55 


5.43 


0.600 


5292. 


-0.676 


4.73 


6.33 


6.61 


6.17 


5.76 


5.20 


4.84 


4.77 


0.700 


5914. 


-0.265 


4.58 


5.30 


5.51 


5.17 


4.83 


4.40 


4.10 


4.07 


0.750 


6282. 


-0.015 


4.47 


4.68 


4.87 


4.57 


4.28 


3.90 


3.64 


3.63 


0.800 


6688. 


0.334 


4.25 


3.81 


3.96 


3.72 


3.47 


3.17 


2.95 


2.96 



Table III. Same as in Table 2 for [M/H] = -1.5 



m 


Teff 


log L 


log g 


Mbol 


Mv 


Mr 


Mi 


Mj 


Mh 


Mk 


0.083 


2194. 


-3.844 


5.51 


14.25 


19.42 


16.32 


13.87 


12.06 


12.24 


12.42 


0.085 


2519. 


-3.559 


5.48 


13.54 


17.50 


15.20 


13.05 


11.43 


11.38 


11.53 


0.090 


2938. 


-3.207 


5.42 


12.66 


15.21 


13.78 


12.13 


10.73 


10.37 


10.37 


0.100 


3255. 


-2.910 


5.34 


11.92 


13.71 


12.61 


11.37 


10.13 


9.63 


9.52 


0.110 


3405. 


-2.745 


5.30 


11.50 


13.04 


12.02 


10.94 


9.77 


9.25 


9.10 


0.130 


3586. 


-2.523 


5.24 


10.95 


12.25 


11.29 


10.38 


9.28 


8.75 


8.57 


0.150 


3691. 


-2.366 


5.20 


10.56 


11.75 


10.81 


9.98 


8.93 


8.40 


8.21 


0.200 


3845. 


-2.082 


5.11 


9.84 


10.89 


10.00 


9.25 


8.26 


7.74 


7.55 


0.300 


4008. 


-1.724 


5.00 


8.95 


9.85 


9.03 


8.34 


7.42 


6.91 


6.72 


0.400 


4134. 


-1.471 


4.92 


8.32 


9.11 


8.35 


7.70 


6.82 


6.31 


6.14 


0.500 


4462. 


-1.116 


4.80 


7.43 


8.00 


7.36 


6.80 


6.01 


5.51 


5.39 


0.600 


5167. 


-0.702 


4.72 


6.40 


6.66 


6.21 


5.79 


5.20 


4.81 


4.75 


0.700 


5792. 


-0.296 


4.58 


5.38 


5.57 


5.21 


4.87 


4.42 


4.11 


4.08 


0.750 


6132. 


-0.059 


4.47 


4.79 


4.96 


4.65 


4.34 


3.95 


3.67 


3.66 


0.800 


6483. 


0.266 


4.27 


3.98 


4.12 


3.86 


3.59 


3.26 


3.02 


3.02 
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Table IV. Sam 


i.c^ us n.i. . 






1 — 












m 


Teff 


log L 


loK K 


Mbol 


Mv 


Mr 


Mi 


Mj 


Mh 


Mk 


u.uoo 




-3.751 


5.49 


14.02 


19.24 


16.30 


13.83 


11.81 


11.82 


11.95 


u.uoo 




-3.524 


5.46 


13.45 


17.61 


15.32 


13.12 


11.34 


11.15 


11.25 


u.uyu 




-3.214 


5.40 


12.68 


15.47 


13.98 


12.23 


10.74 


10.31 


10.27 


U.iUU 


ozUU. 


-2.931 


5.34 


11.97 


13.93 


12.79 


11.46 


10.16 


9.63 


9.49 


0.110 


3343. 


-2.769 


5.29 


11.56 


13.24 


12.20 


11.03 


9.81 


9.26 


9.09 


0.130 


3512. 


-2.552 


5.23 


11.02 


12.44 


11.46 


10.47 


9.33 


8.77 


8.58 


0.150 


3624. 


-2.392 


5.19 


10.62 


11.90 


10.95 


10.05 


8.97 


8.42 


8.22 


0.200 


3780. 


-2.107 


5.10 


9.91 


11.02 


10.11 


9.32 


8.30 


7.77 


7.57 


0.300 


3945. 


-1.747 


4.99 


9.01 


9.97 


9.12 


8.41 


7.46 


6.94 


6.74 


0.400 


4067. 


-1.492 


4.92 


8.37 


9.22 


8.42 


7.76 


6.85 


6.33 


6.16 


0.500 


4377. 


-1.139 


4.79 


7.49 


8.12 


7.44 


6.86 


6.04 


5.51 


5.39 


0.600 


5068. 


-0.727 


4.71 


6.46 


6.75 


6.27 


5.84 


5.22 


4.81 


4.74 


0.700 


5686. 


-0.326 


4.58 


5.46 


5.64 


5.27 


4.92 


4.44 


4.12 


4.09 


0.750 


6006. 


-0.099 


4.47 


4.89 


5.05 


4.72 


4.41 


3.99 


3.71 


3.69 


0.800 


6315. 


0.199 


4.29 


4.14 


4.29 


4.00 


3.72 


3.36 


3.11 


3.10 



Table V. Same as in Table 2 for [M/H] = -1.0 



m 


Teff 


log L 


log g 


Mbol 


Mv 


Mr 


Mi 


Mj 


Mh 


Mk 


0.083 


2359. 


-3.660 


5.45 


13.79 


19.10 


16.42 


13.93 


11.60 


11.35 


11.42 


0.085 


2550. 


-3.492 


5.43 


13.37 


17.81 


15.58 


13.31 


11.26 


10.89 


10.91 


0.090 


2840. 


-3.231 


5.38 


12.72 


15.94 


14.35 


12.44 


10.75 


10.24 


10.15 


0.100 


3107. 


-2.967 


5.32 


12.06 


14.35 


13.15 


11.64 


10.21 


9.63 


9.46 


0.110 


3245. 


-2.808 


5.28 


11.66 


13.61 


12.52 


11.19 


9.87 


9.28 


9.09 


0.130 


3406. 


-2.593 


5.22 


11.12 


12.74 


11.74 


10.61 


9.39 


8.81 


8.60 


0.150 


3505. 


-2.439 


5.18 


10.74 


12.20 


11.23 


10.20 


9.04 


8.46 


8.25 


0.200 


3674. 


-2.150 


5.10 


10.02 


11.25 


10.33 


9.45 


8.37 


7.82 


7.60 


0.300 


3841. 


-1.787 


4.99 


9.11 


10.17 


9.28 


8.51 


7.52 


6.98 


6.78 


0.400 


3959. 


-1.530 


4.91 


8.46 


9.42 


8.57 


7.86 


6.91 


6.38 


6.18 


0.500 


4235. 


-1.185 


4.78 


7.60 


8.35 


7.60 


6.98 


6.11 


5.55 


5.42 


0.600 


4867. 


-0.782 


4.69 


6.60 


6.96 


6.42 


5.97 


5.28 


4.81 


4.74 


0.700 


5478. 


-0.392 


4.58 


5.62 


5.81 


5.41 


5.04 


4.52 


4.17 


4.13 


0.750 


5776. 


-0.181 


4.49 


5.09 


5.25 


4.89 


4.56 


4.10 


3.79 


3.77 


0.800 


6055. 


0.077 


4.34 


4.45 


4.59 


4.27 


3.96 


3.56 


3.28 


3.27 



